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Abstract 

Radiatively-driven flow in a luminous disk is examined in the subrelativistic regime of (v/c) 1 , taking 
account of radiation transfer. The flow is assumed to be vertical, and the gravity and gas pressure are 
ignored. When internal heating is dropped, for a given optical depth and radiation pressure at the flow 
base (disk "inside"), where the flow speed is zero, the flow is analytically solved under the appropriate 
boundary condition at the flow top (disk "surface"), where the optical depth is zero. The loaded mass 
and terminal speed of the flow are both determined by the initial conditions; the mass-loss rate increases 
as the initial radiation pressure increases, while the flow terminal speed increases as the initial radiation 
pressure and the loaded mass decrease. In particular, when heating is ignored, the radiative flux F is 
constant, and the radiation pressure Pq at the flow base with optical depth tq is bound in the range of 
2/3 < cP /F < 2/3 + tq. In this case, in the limit of cPq/F = 2/3 + To, the loaded mass diverges and the 
flow terminal speed becomes zero, while, in the limit of cPq/F = 2/3, the loaded mass becomes zero and 
the terminal speed approaches (3/8)c. which is the terminal speed above the luminous flat disk under an 
approximation of the order of (v/c) 1 . We also examine the case where heating exists, and find that the 
flow properties are qualitatively similar to the case without heating. 
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— relativity — X-rays: stars 



1. Introduction 

Accretion disks are tremendous energy sources in the 
active universe: in young stellar objects, in cataclysmic 
variables, in galaxtic X-ray sources and microquasars, and 
in active galaxies and quasars (see Kato et al. 1998 for 
a review). In particular, in a supercritical accretion disk, 
the mass-accretion rate highly exceeds the critical rate, 
the disk local luminosity exceeds the Eddington one, and 
the mass loss from the disk surface driven by radiation 
pressure takes place. 

Such a radiatively driven outflow from a luminous ac- 
cretion disk has been extensively studied in the con- 
text of models for astrophysical jets by many researchers 
(Bisnovatyi-Kogan, Blinnikov 1977; Katz 1980; Icke 1980; 
Melia, Konigl 1989; Misra, Melia 1993; Tajima, Fukue 
1996, 1998; Watarai, Fukue 1999; Hirai, Fukue 2001; 
Fukue ct al. 2001; Orihara, Fukue 2003), as on-axis jets 
(Icke 1989; Sikora et al. 1996; Renaud, Henri 1998; Luo, 
Prothcroe 1999; Fukue 2005a), as outflows confined by a 
gaseous torus (Lyndcn-Bcll 1978; Davidson, McCray 1980; 
Sikora, Wilson 1981; Fukue 1982), or as jets confined by 
the outer flow or corona (Sol et al. 1989; Marcowith et al. 
1995; Fukue 1999), and as numerical simulations (Eggum 
et al. 1985, 1988). In almost all of these studies, however, 
the disk radiation field was treated as an external field, 
and the radiation transfer was not solved. 

The radiation transfer in the disk, on the other hand, 
has been investigated in relation to the structure of a 
static disk atmosphere and the spectral energy density 



from the disk surface (e.g., Meyer, Meyer-Hofmeister 1982; 
Cannizzo, Wheeler 1984; Shaviv, Wehrse 1986; Adam et 
al. 1988; Hubeny 1990; Ross et al. 1992; Artemova et 
al. 1996; Hubeny, Hubeny 1997, 1998; Hubeny et al. 
2000, 2001; Davis et al. 2005; Hui et al. 2005; see 
also Mincshige, Wood 1990). In these studies, however, 
the vertical movement and mass loss were not considered. 
Moreover, the relativistic effect for radiation transfer was 
not considered in the current studies. 

In this paper we first examine the radiatively driven ver- 
tical flow - moving photosphere - in a luminous flat disk 
within the framework of radiation transfer in the subrel- 
ativistic regime of (v/c) 1 , where we incorporate the effect 
of radiation drag up to the first order of the flow velocity. 
At the first stage, we ignore the gravity of the central ob- 
ject as well as the gas pressure, although internal heating 
is considered in a limited way. We then consider analyt- 
ical solutions for such a radiative flow, and examine the 
flow properties. 

In the next section we describe basic equations in the 
vertical direction. In section 3 we examine the radia- 
tive flow without internal heating, while the radiative flow 
with internal heating is discussed in section 4. The final 
section is devoted to concluding remarks. 

2. Basic Equations 

Let us suppose a luminous flat disk, deep inside which 
gravitational or nuclear energy is released via viscous heat- 
ing or other processes. The radiation energy is trans- 
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ported in the vertical direction, and the disk gas, itself, 
also moves in the vertical direction due to the action of ra- 
diation pressure (i.e., plane-parallel approximation). For 
simplicity, in the present paper, the radiation field is con- 
sidered to be sufficiently intense that both the gravita- 
tional field of, e.g., the central object and the gas pressure 
can be ignored: tenuous cold normal plasmas in the super- 
Eddington disk, cold pair plasmas in the sub-Eddington 
disk, or dusty plasmas in the sub-Eddington disk. As 
for the order of the flow velocity v, we consider the sub- 
relativistic regime, where the terms of the first order of 
(v/c) arc retained, in order to take account of radiation 
drag. Under these assumptions, the radiation hydrody- 
namic equations for steady vertical (z) flows are described 
as follows (Kato et al. 1998). 
The continuity equation is 

pv = J (= const.), (1) 

where p is the gas density, v the vertical velocity, and J 
the mass-loss rate per unit area. The equation of motion 
is 

dv K a bs + K sca 



[F-(E + P)v], 



(2) 



dz c 

where K a b s and « sca are the absorption and scattering 
opacities (gray), E the radiation energy density, F the 
radiative flux, and P the radiation pressure. In a gas- 
pressureless approximation, the energy equation is re- 
duced to 



= q + - P [j - CK a b s S + K 



2Fv 



abs " 



(3) 



where q + is the heating and j is the emissivity. For radi- 
ation fields, we have 



dF 

dz 



= P 



Fv 

J CK absF + (K a bs — K sca ) 

C 



(4) 



^ = BZ (j - cn ahs E) / abs + Ksca [F — (E + P)v],(5) 



and the closure relation, 
1 4Fv 



(6) 



Eliminating j and E with the help of equations (3) and 
(6), and introducing the optical depth by 

dr = -(K abs + K sca )pdz, 

we can rearrange the basic equations up to the order of 

{v/c) 1 as 



dv 

cJ— =-(F- APv) , 



(^abs H~ ^scel)J 



dF 
~ch r 
dP 

~d7 

dz 



q 



( K abs + ^sca)P 

F - 4Pv, 



F V ~, 
c 



(8) 
(9) 
(10) 
(11) 



Finally, integrating the sum of equations (8) and (10) 
gives the momentum conservation in the present approx- 
imation, 



Jv + P = K (= const.), 



(12) 



and similarly from equations (8) and (9) we have the en- 
ergy conservation, 



1 



~-Jv 2 + F- 



-dr, 



(13) 



(^abs + Ksca)/? 

where the first term on the left-hand side is eventually 
dropped, although we retain it here to clarify the physical 
meanings. 

We solved equations (8), (12), and (13) for appropriate 
boundary conditions, and we obtained analytic solutions 
for several special cases. 

As for the boundary conditions, wc imposed the fol- 
lowing cases. At the flow base (disk "inside") with an 
arbitrary optical depth To, the flow velocity is zero and 
the radiation pressure is Pq, where subscript denotes 
the values at the flow base. At the flow top (disk "sur- 
face"), where the optical depth is zero, the radiation fields 
should satisfy the values above the luminous infinite disk. 
Namely, just above the disk with surface intensity Iq, the 
radiation energy density E s , the radiative flux F s , and 
the radiation pressure P s are E s = {2/c)txIq, F s = txIq, 
and P s = (2/3c)7rJ , respectively, where subscript s de- 
notes the values at the flow top. Rigorously speaking, this 
latter boundary condition at the flow top is for a static 
photosphere, and should not be applied to a moving pho- 
tosphere, although we approximately used this condition 
in the present subrclativistic regime. The exact boundary 
conditions are derived and discussed in a separate paper 
for a fully special relativistic case. 

As shown explicitly later, in order for the flow to exist, 
the radiation pressure Po at the flow base should be re- 
stricted in some ranges. The lower limit is set due to the 
negative gradient of the radiation pressure toward the flow 
top. The upper limit is set when the loaded mass diverges. 
Beyond the upper limit, the hydrostatic balance without 
mass loss may be established in the vertical direction. 

3. Radiative Flow Without Heating 

We first examine the case, where there is no viscous 
(7) or nuclear heating: q + = 0. In the luminous accretion 



disk, it is usually supposed that the heating processes take 
place deep inside the disk. Hence, the assumption without 
heating approximately describes the radiative flow in or 
near the surface envelope of the disk. 

In this case and under the present subrelativistic 
regime, equation (13) means that the radiative flux F is 
conserved along the flow, 

F = F S . (14) 

Using the boundary conditions at the flow base (v = 0), 
equation (12) is expressed as 



Jv + P = P 



(15) 
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Hence, equation (8) becomes 
dv 



c.J 



{F s - AP {) v) . 



(16) 



This equation can be analytically solved to yield the so- 
lution for v, 



AP 



1 



(17) 



Thus, the radiative flow from the luminous disk without 
heating is expressed in terms of the boundary values and 
the mass-loss rate. In addition, the flow velocity v s at the 
flow top (r = 0) is 

I-e~^r T °Y (18) 




Using the boundary condition at the flow top, we fur- 
ther impose a condition on the values at boundaries. As 
already mentioned, at the flow top cP = cP s = (2/3)F s . 
Hence, inserting boundary values into momentum equa- 
tion (15), we have the following relation: 



J- 



4P 



_4P 
1-e cJ 



TO 



3c 



(19) 



That is, for given To and Pq at the flow base, the mass-loss 
rate J is determined in units of F s /c 2 , as an eigen value. 

In figure 1 we show the mass-loss rate J (solid curves) 
and the terminal velocity v s (dashed curves) at the flow 
top as a function of Po for several values of To at the flow 
base. The quantities are normalized in units of c and F s . 

As can be seen in figure 1, the mass- loss rate increases 
as the initial radiation pressure increases, while the flow 
terminal speed increases as the initial radiation pressure 
and the loaded mass decrease. 

Moreover, as can be seen in figure 1, and easily shown 
from equation (19), in order for the flow to exist, the ra- 
diation pressure Pq at the flow base is restricted in some 
range, 



2 cP 2 

3 < ~F < 3 



•70- 



(20) 



At the upper limit of cPq/F s = 2/3 + t , the loaded mass 
diverges and the flow terminal speed becomes zero. On 
the other hand, at the lower limit of cP /F s = 2/3, the 
pressure gradient vanishes, the loaded mass becomes zero 
and the terminal speed approaches (3/8)c. This is just the 
terminal speed of the gas, driven by the radiation pressure 
under the radiation drag, above the luminous Hat infinite 
disk in the present approximation of the order of (v/c) 1 
(cf. Icke 1989). 

In figure 2 we show the flow velocity v (solid curves) and 
the radiation pressure P (dashed curves) as a function of 
the optical depth t for several values of Pq at the flow 
base in a few cases of To. The quantities are normalized 
in units of c and F s . 

When the initial radiation pressure Pq at the flow base 
is large, the pressure gradient between the flow base and 
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Fig. 1. Mass-loss rate J (solid curves) and terminal velocity 
v B (dashed curves) at the flow top as a function of Pq for 
several values of To at the flow base: (a) to = 1 and (b) tq = 10. 
The quantities are normalized in units of c and _F S . 

the flow top is also large. As a result, the loaded mass 
J also becomes large, but the flow final speed v s is small 
due to momentum conservation (12). When the initial 
radiation pressure Pq [> (2/3c)F s ] is small, on the other 
hand, the pressure gradient becomes small, and the loaded 
mass is also small, but the flow final speed becomes large. 
In the limiting case of Pq = (2/3c)F s , the loaded mass 
J becomes zero, but the terminal speed v s is (3/8)c, as 
already mentioned. 

Finally, inserting equation(17) into equation (11), and 
integrating the resultant equation, we obtain the height z 
as a function of the optical depth t: 



(^abs "I" Ksca)J 



F s cJ 
^P \ iT °- T) -4p- 



(21) 



where we use the boundary condition at the flow base: 
z = at t = tq. In addition, the height z s of the flow top 
(disk "surface" ) is also expressed as 



(« 



abs 



4cP n 



To + 3 



cPo 
F* 



(22) 
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Fig. 2. Flow velocity v (solid curves) and radiation pressure 
P (dashed curves) as a function of the optical depth r for 
several values of Po at the flow base in a few cases of to: (a) 
to = 1 and (b) to = 10. The quantities are normalized in units 
of c and F B . From top to bottom of v and from bottom to 
top of P, the values of P are 0.67, 0.8, 1, 1.2, 1.4 in (a), and 
0.67, 1, 2, 5 in (b). 

where we use equation (18). Hence, at the upper limit of 
cPq/F s = 2/3 + t , where the loaded mass diverges, the 
height of the flow/disk becomes zero. On the other hand, 
at the lower limit of cPq/ F s = 2/3, where the loaded mass 
becomes zero, the height of the flow/disk becomes infinite. 

It should be noted that, at the upper limit of cPo/F s = 
2/3 + tq with a small terminal speed, 



T 



err 



3 +T ° 



(23) 



where T c is the temperature at the flow base (i.e., Po = 
aT^/3), and T e g the effective temperature at the flow top 
(i.e., F s = crT^ s ). This is the Milne approximation. 



4. Radiative Flow With Heating 

We next examine the effect of heating: q + ^ 0. That 
is, we consider the radiative flow from deep inside of the 
disk. 

If the quantity q + / p is a function of r, the heating term 
can be integrated [see equations (9) or (13)]. In order 



for the integration to be performed analytically, several 
functional forms need to be assumed. One simple case, 
which we assume in the present paper, is a following form: 

* = q (= const.). (24) 

This is quite reasonable since the heating q + may be pro- 
portional to the density p. 

In this case, under the present subrelativistic regime 
and using the boundary conditions at the flow base, equa- 
tions (8), (13), and (12) become, respectively, 



cJ^ = -(F-4P Q v) 1 



F = <?o(to - r) 
P = Pa- Jv, 



— [TQ - T) , 
TO 



(25) 

(26) 
(27) 



where F s (= go r o) is again the radiative flux at the flow top 
(disk "surface"). Here, we impose the boundary condition 
on F at the flow base as F = at r = ro. 

The above equation (25) can be analytically solved to 
yield the solution for v: 



F 



4P ro 



(ro-r)- 



c,JF s 
(4Po) 2 r 



1-e cJ 



4i V-T ) 



(28) 



Thus, the radiative flow from the luminous disk with sim- 
ple heating is also expressed in terms of the boundary val- 
ues and the mass-loss rate. In addition, the flow velocity 
v s at the flow top (t = 0) is 

4P . 



4P 



cJF s 
(4Po) 2 T 



1 



cJ 



T 



(29) 



Using the boundary condition at the flow top, we fur- 
ther impose a condition on the values at the boundaries. 
In this case with simple heating, we have the following 
relation: 

4Pn, 
1-e cJ 



J 



cJ*R I -4P, 



T 



4P (4P ) 2 r 



3c 



-Fs = Pn 



(30) 



That is, for given To and Po at the flow base, the mass-loss 
rate J is determined in units of F s /c 2 . 

In figure 3 we show the mass-loss rate J (solid curves) 
and the terminal velocity v s (dashed curves) at the flow 
top as a function of Po for several values of To at the flow 
base in the case with simple heating. The quantities are 
normalized in units of c and F s . 

The properties are qualitatively same as those without 
heating. The mass-loss rate increases as the initial ra- 
diation pressure increases, while the flow terminal speed 
increases as the initial radiation pressure and the loaded 
mass decreases. Moreover, in order for flow to exist, the 
radiation pressure Po at the flow base is restricted in some 
range: 

2 cP 2 1 

3 < ^ < 3 + 2 r °- (31) 




Fig. 3. Mass-loss rate J (solid curves) and terminal velocity 
v B (dashed curves) in the case with simple heating: (a) To = 1 
and (b) to = 10. The quantities are normalized in units of c 
and F s . 



At the upper limit, the loaded mass diverges and the flow 
terminal speed becomes zero. On the other hand, at the 
lower limit of cPq/F s = 2/3, the loaded mass becomes zero 
and the terminal speed approaches (3/8)c. 

In figure 4 we show the flow velocity v (solid curves) and 
the radiation pressure P (dashed curves) as a function of 
the optical depth r for several values of Pq at the flow 
base in a few cases of To in the case with simple heating. 
The quantities are normalized in units of c and F s . 

In this case with simple heating, the velocity distribu- 
tion is somewhat different from that in the case without 
heating. In the case without heating, the flow is rather 
quickly accelerated to be a terminal speed, since the ra- 
diative flux is constant. In the case with heating, on the 
other hand, the flow is gradually accelerated, since the 
radiative flux is increasing from the flow base toward the 
flow top due to heating. 

Finally, inserting equation(17) into equation (11), and 
integrating the resultant equation, we obtain the height z 
as a function of the optical depth t: 



(^abs ~t~ K SCSL )J 



Fig. 4. Flow velocity v (solid curves) and radiation pressure 
P (dashed curves) in the case with simple heating: (a) to = 1 
and (b) tq = 10. The quantities are normalized in units of c 
and F B . From top to bottom of v and from bottom to top of 
P, the values of P are 0.7, 0.8, 0.9, 1, 1.1 in (a), and 0.7, 1, 
2, 5 in (b). 



4cP 



cJ 

1F Q 



2r 
" 1 

TQ 



cJ 
4P r 



(ro-r) 



cJ 



(r-T ) 



(32) 



where we use the boundary condition at the flow base: 
z = at t = tq. In addition, the height z s of the flow top 
(disk "surface" ) is also expressed as 



(«abs ' 



0-7 



4cP n 



cPo 



(33) 



where we use equation (29). Hence, at the upper limit of 
cPo/F s = 2/3 + t /2, where the loaded mass diverges, the 
height of the flow/disk becomes zero. On the other hand, 
at the lower limit of cPq/F s = 2/3, where the loaded mass 
becomes zero, the height of the flow/disk becomes infinite. 

It should be noted that, at the upper limit of cPq/F s = 
2/3 + tq/2 with small terminal speed, 



(34) 
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where T c is the temperature at the flow base (i.e., Po = 
aT^/3), and T e ff the effective temperature at the flow top 
(i.e., F s = crT c 4 ff ). This is just the modified Milne approx- 
imation. 

5. Concluding Remarks 

In this paper we have examined the radiative flow from a 
luminous disk, while taking account of the radiative trans- 
fer, in the subrelativistic regime of (v/c) 1 . The flow ve- 
locity, the radiation pressure distribution, and other quan- 
tities are analytically solved as a function of the optical 
depth for the cases without and with heating. At the flow 
base (disk "inside"), where the flow speed is zero, the ini- 
tial optical depth is To and the initial radiation pressure 
is Po; in the usual accretion disk these quantities are de- 
termined in terms of the central mass, the mass-accretion 
rate, and the viscous process as a function of the radius. 
That is, the optical depth is determined by the disk sur- 
face density, the radiative flux is given by the central mass 
and the accretion rate, and the radiation pressure relates 
to the disk internal structure. At the flow top (disk "sur- 
face"), where the optical depth is zero, the radiation fields 
are assumed to be those above a static flat disk with uni- 
form intensity. In order to match this boundary condition, 
one relation is imposed on the boundary quantities, and 
the mass-loss rate J and terminal speed v s are both deter- 
mined by the initial conditions To and Po, as eigenvalues. 
In particular, in order for the flow to exist, the initial radi- 
ation pressure is bound in some range: at the upper limit, 
the loaded mass diverges and the flow terminal speed be- 
comes zero, while, at the lower limit, the loaded mass 
becomes zero and the terminal speed approaches (3/8)c, 
which is the terminal speed above the luminous disk under 
the approximation of the order of (v/c) 1 . 

Among various types of astrophysical jets and outflows, 
some are believed to be accelerated by continuum radi- 
ation, and to which the present mechanism can be ap- 
plied. For example, in supersoft X-ray sources, including 
accreting white dwarfs, high-velocity mass outflows were 
reported (e.g., Cowley et al. 1998). The wind velocity 
measured in these objects is about 3000-5000 km s" 1 , 
which is of the order of the escape velocity of white dwarfs. 
In addition, the line profile depends on the inclination an- 
gle, and the P Cyg profiles are seen. Hence, in supersoft 
X-ray sources, the outflows should blow off from the inner- 
most accretion disk. However, the gas temperature of the 
inner disk is so high that the line-driven mechanism can- 
not operate, and therefore, mass outflows in these objects 
would be driven by continuum radiation. 

In SS 433, a prototype of astrophysical jets, the mass- 
accretion rate highly exceeds the Eddington one, and twin 
jets with 0.26 c would be accelerated by radiation pres- 
sure (see, e.g., Cherepashchuk et al. 2005 for a recent ob- 
sevation). The luminosity of microquasar GRS 1915+105 
exceeds the Eddington one (^Edd) during its outburst. In 
quiscent phase, its luminosity never drop below ~ 0.3LEdd 
(Done, Gierlihski 2004). 

In broad absorption- line quasars (BAL QSOs), strong 



broad absorption lines in the UV resonance are always 
blucshiftd relative to the emission-line rest frame, which 
indicates the presence of outflows from the nucleus, with 
velocities as large as 0.2 c. These moderately high- velocity 
outflows are supposed to be accelerated by line or contin- 
uum radiation (e.g., Murray et al. 1995; Proga et al. 2000; 
Proga, Kallman 2004). In these objects, the outflows must 
consist of dense clouds, or be shielded from the central 
strong source, in order to operate a line-driven mecha- 
nism. Otherwise, the material is highly ionized by the 
central source, and, instead of line-acceleration, contin- 
uum acceleration may work (see, e.g., Chclouche, Netzer 
2003). 

Finally, narrow-line Seyfert I galaxies and luminous 
quasars are also believed to harbor supercritically accret- 
ing black holes, and high- velocity outflows are reported as 
warm absorbers (cf. Blustin et al. 2005). In these objects, 
radiative acceleration by continuum may play a dominant 
role. 

The radiative flow investigated in the present paper 
must be quite fundamental problems for accretion disk 
physics and astrophysical jet formation, although the 
present paper is only the first step, and there are many 
simplifications at the present stage. 

For example, we have ignored the gravitational field 
produced by the central object. This means that the flow 
considered in the present paper would correspond to nor- 
mal plasmas in the super-Eddington disk, pair plasmas in 
the sub-Eddington disk, or dusty plasmas in the luminous 
disk. In other cases, or even for normal plasmas in the 
super-Eddington disk, the gravitational field would affect 
the flow properties. In particular, the gravitational field 
in the vertical direction is somewhat complicated (e.g., 
Fukuc 2002, 2004), and the influence of the gravitational 
field is important. The effect of gravity will be considered 
in a seperate paper (Fukue 2005b). 

We also ignored the gas pressure. This means that the 
radiation field is sufficiently intense. In general cases, 
where the gas pressure is considered, there usually ap- 
pear sonic points (e.g., Fukue 2002, 2004), and the flow 
is accelerated from subsonic to supersonic. The cross sec- 
tion of the flow also has a similar influence. In this paper 
we consider a purely vertical flow, and the cross section 
of the flow is constant. If the cross section of the flow 
increases along the flow, the flow properties, such as the 
transonic nature, would be influenced. 

In the present paper, we have considered the internal 
heating processes in a simple form. Since the heating 
prcesses should couple with the viscous process or other 
heating processes, such as a nuclear process, we should 
treat internal heating more carefully. 

In order to take account of the effect of radiation drag, 
we consider the subrelativistic regime, where the terms 
of the first order of (v/c) are retained. As a result, the 
terminal speed (3/8)c for the subrelativistic regime ap- 
pears in the extreme case. This is quite consistent with 
the previous knowledge. As is well known, the terminal 
speed above the luminous flat disk in the full relativity 
is [(4-\/7)/3]c ~ 0.45c (Icke 1989). To consider highly 
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relativistic cases, all of the terms up to {v/c) 2 should be 
retained (cf. Kato et al. 1998). A related problem is the 
boundary condition at the flow top. In the present paper, 
we impose the boundary condition such that the radia- 
tion fields at the flow top are those above the luminous 
flat disk. Rigorously speaking, this boundary condition 
should be modified, since the disk gas, itself, is now mov- 
ing. In the highly relativistic case, the deviation would be 
serious. The exact boundary conditions are derived and 
discussed in a separate paper for a fully special relativistic 
case (Fukue 2005c). 

This work has been supported in part by a Grant-in- Aid 
for the Scientific Research (15540235 J.F.) of the Ministry 
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